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I. INTRODUCTION 

The experimental demonstration of Bose-Einstein condensation in atomic vapors ^-^ 



Q I has rapidly lead to spectacular new advances in atom optics. In particular, it has enabled 



its extension from the linear to the nonlinear regime, very much like the laser lead to the 

development of nonlinear optics in the 1960s. It is now well established that two-body 

^ ■ collisions play for matter waves a role analogous to that of a Kerr nonlinear crystal in 

^ ! optics. In particular, it is known that the nonlinear Schrodinger equation which describes 

^ I the condensate in the Hartree approximation supports soliton solutions. For the case of 

^ I repulsive interactions normally encountered in BEG experiments, the simplest solutions are 

On ' dark solitons, that is, 'dips' in the density profile of the condensate. These dark solitons have 

O 



been recently demonstrated in two experiments [^,0 which appear to be in good agreement 



O . with the predictions of the Gross-Pitaevskii equation. 



While very interesting from a fundamental physics point-of-view, dark solitons would 

appear to be of limited potential for applications such as atom interferometry, where it is 

desirable to achieve the dispersionless transport of a spatially localized ensemble of atoms, 

^ ■ rather than a 'hole'. In that case, bright solitons are much more interesting. However, the 

O . problem is that large condensates are necessarily associated with repulsive interactions, for 

which bright solitons might seem impossible since the nonlinearity cannot compensate for 

the kinetic energy part (diffraction) in the atomic dynamics. While this is true for atoms in 

/\ • free space, this is however not the case for atoms in suitable potentials, eg. optical lattices. 

c^ ■ This is because in that case, it is possible to tailor the dispersion relation of the atoms in 

such a way that their effective mass becomes negative. For such negative masses, a repulsive 

interaction is precisely what is required to achieve soliton solutions. This result is known 

from nonlinear optics, where such soliton solutions, called gap solitons, have been predicted 

and demonstrated [§-0. 

The article is organized as follows: In Section ^ we briefly review the basic formalism 
for nonlinear atom optics with a focus on Bose condensed systems. This leads to the Gross- 
Pitaevskii equation, the nonlinear mean-field equation of motion for the atomic condensate. 
In Section PTI we introduce the bright and dark soliton solutions of these equations and 



motivate the concept of gap solitons. In Section |^ we then describe a specific model for 
gap solitons in a spinor Bose-Einstein condensate. We demonstrate how to launch and 
control these solitons in Section 0. Finally, Section ^ is a summary and conclusion. 



II. NONLINEAR ATOM OPTICS 



Early experiments and theories concerning atom optics considered the low-density regime 
where atom-atom interactions are negligible. However, as the prospect of atomic BEC be- 
came more tangible prior to 1995, Lenz et al. [[10| , pT| and Zhang et al. [|1^ theoretically 



explored the high density regime where atom-atom interactions become relevant. In partic- 
ular, they view the cold atomic collisions as nonlinear mixing processes analogous to those 
in nonlinear optics, and the new area of nonlinear atom optics was born. These initial 
theoretical works considered the dipole-dipole interaction between atoms, and derived a 
mean-field equation for the macroscopic wave function or order parameter for the system 
of cold atoms analogous to the well-known nonlinear Schrodinger equation from nonlinear 
optics |13|, or the Gross-Pitaevskii equation (GPE) from BEC theory |14,15|. It was shown 



that these equations admit soliton solutions under suitable conditions |T^,[T2|. The experi- 
mental demonstrations of atomic BEC in 1995 in both Rubidium [jl| and Sodium [Q ushered 
in the era of experimental nonlinear atom optics, where the cold atomic collisions are now 
mediated via the Van der Walls interaction, with experimental demonstrations of atomic 
four- wave mixing [16|, dark atomic solitons 0,^, atomic vortices [|l^,|l^, and mode-locking 
of an atom laser 



19 



To set the stage for our discussion, we first briefly review the physics and derivation of 
the mean-field GPE for an atomic BEC, which underlies nonlinear atom optics. We assume 
that an ensemble of ultracold atoms of mass m experience a single-particle Hamiltonian of 
the form 



Ho = -I^V^ + V{K) 
2m 



(1) 



where the external potential ^(R) could be for example an optical dipole potential pO|. In 
the language of second-quantization this translates to the Hamiltonian operator 



no= |ciV^^(R)ifo^(R) 



(2) 



where \i/^(R) and ^I/(R) denote the creation and annihilation operators (in the Schrodinger 
picture) for an atom with center-of-mass position vector R. Assuming bosonic atoms, the 
field operators obey the commutation relation [^(R), \&t(R')] = 5(R — R'). Next we add 
atom-atom interactions by introducing a two-body potential V(R, R'), which in second- 
quantization leads to the total Hamiltonian operator Ti. = T^o + V with 



V 



d^rd^r'^\R)^\K')V(R, R0^(R')^(R) • 



(3) 



For a system of N cold atoms, the many-body state vector may be related to the A^- 
particle Schrodinger wave function using the general relation 



Mt)) = ^ / ^'^1 • • • d^r^M^u • • • , Rtv, t)^^(R7v) . . . ^^(Ri)|0) . 



(4) 



Since for a BEC at zero temperature the atoms are predominantly in the same Bose- 
condensed state, we employ the Hartree approximation in which the A-particle wave function 
{/)Ar(Ri, . . . , Rjv, t) is written in the form of a product 



N 



'JN 



;Ri,...,Rjv,t) = n^(R^.^) 



(5) 



j=i 



where A/iV'i9(R, t) = (0Ar_i|^(R)|0Ar) is the effective single-particle, or Hartree, wave func- 
tion, and is assumed normalized [^. We note that the Hartree approximation automati- 
cally satisfies the symmetry condition for bosons that interchanging any two particle indices 
produces the same wave function. To obtain an equation of motion for the effective single- 
particle wave function, we employ Dirac's variational principle [B^] 



(5^*(R,t) 



VI.4 



-Ho + V 



JN 







(6) 



which yields the GPE 



on f 

ih— = Hod + {N -I) dVV(R, R') |i?(R', t) |279(R, t) . 



(7) 



Here the nonlinear term describes the mean- field effect of the (A^ — 1) other atoms on the 
effective single-particle evolution. 

For the particular case of atomic BECs the cold atomic collisions may be treated in the 
s-wave scattering limit. The effective two-body potential can then be approximated by a 



contact potential of the form [23,24 



y (R, R') = f/o(5(R - R') 



(8) 



where the coefficient Uq = AtxU asc/m^ with asc the s-wave scattering length. The GPE 
equation (1^) for an atomic BEG then becomes 



ifi 



'dt 



2m 



V¥ + 1/(R)^ + A^[/o|^|¥, 



(9) 



which applies for large A^. This equation forms the basis of the following discussion of atomic 
solitons. Here we have concentrated on the case of a single atomic species, but the formalism 



is straightforward to generalize to multi-component or vector condensates p5| , |26[| involving, 
for example, several different Zeeman sublevels of a given atom. 



III. ONE-DIMENSIONAL ATOMIC SOLITONS 



A. One-dimensional GPE 



In this Section we briefly discuss the one-dimensional (ID) soliton solutions of the GPE 
(§). In particular, we assume strong transverse confinement of the BEG in the {X, Y) plane 
such that the transverse mode Ug{X, Y) is the ground state of the potential V{X, Y), 






Ug + V{X,Y)U, 



(10) 



This may be realized, for example, using a red-detuned optical dipole trap around the 
focus of a Gaussian laser beam centered on the origin |^|, thereby giving approximately 
harmonic transverse confinement. Thus, the BEG becomes essentially one-dimensional with 
cross-sectional area At = J dXdY\ug{X,Y)/ug{0,0)\'^ transverse to the unbound Z-axis. 
We remark that the quasi-lD nature of the system does not preclude BEG along with the 
associated GPE description, as discussed by Petrov et al. p^. Setting 



V Ug{0,0) J 



t) 



(11) 



SO that if is normalized to the number of atoms N, and projecting onto the ground transverse 
mode, the reduced ID GPE becomes (we ignore the variation of the optical potential along 
the Z-axis for simplicity) 



ih 



dip 



n" 9V 



dt 2m dZ'^ 

with U = UoA^^ J dXdY\ug{X,Y)/ug{0,0)\\ 



+ U\if\^<f, 



(12) 



B. Dark solitons 



Soliton solutions are non-spreading solutions of Eq. (|I2D which preserve their shape under 
The kinetic energy term in the GPE tends to spread wavepackets, and 



propagation 28 



the nature of the solitons depends on the sign of the nonlinearity. For the case of repulsive 
interactions {ag^ > 0), both the kinetic energy and nonlinearity terms in Eq. (|T2|) tend to 
broaden localized wavepackets, so we do not expect localized, or bright soliton, solutions for 
that case. However, dark solitons describing localized density dips in an otherwise constant 
background can arise and are given by 0,|,^] (with analogous solutions being well known 
in nonlinear fiber optics e.g. 



v{z,t) 



n 



1/2 



\ 



1 2 sech'^ 

^0/ 



{Z - vt) 



1/2N 



j((P(z,t,v)-fj.t/h) 



(13) 



h /mil 



where n is the background density away from the dark soliton core, /i = n|f/|, Zq = 

is the correlation length which determines the width of the soliton core, Vq = J n/m is the 
Bogoliubov speed of sound, v is the dark soliton velocity, whose magnitude is bounded by 
fo, and the soliton phase is given by 



[Z, t,v) = — arctan 




tanh 




(14) 



These solutions show that the dark solitons are characterized by the presence of a phase-step 
S across the localized density dip. It can be related to the velocity v and the density Ubot at 
the bottom of the atomic density dip HJ3]. In particular, one finds the relation 



COS 



V 



Ubot 



n 



(15) 



so that for a stationary soliton v = O,nbot = 0, and there is a. 6 = rr phase-step. Only 
stationary sohtons have a vanishing density at their center, so they are also referred to as 
'black' solitons, whereas for the Ubot > case the expression 'grey' soliton is used. 

Dark matter-wave solitons have recently been observed experimentally in a Na BEG 
by Denschlag et al. [Q, and in a cigar-shaped *^Rb condensate by Burger et al. [^fl. In 
both experiments dark solitons of variable velocity were launched via the phase imprinting 
of a BEG by a light-shift potential. By applying a pulsed laser to only half of the BEG 
and choosing the laser intensity and duration to select a desired phase-step 5, the soliton 
velocity could be selected according to Eq. (^). Fig. |1] shows a schematic of the basic 
phase-imprinting idea for a cigar-shaped BEG, and Fig. ^ shows the experimental results 
of Burger et al. The motion of the density dip is seen for a phase-step of vr and a variety 
of evolution times. Ideally this value of the phase-step should lead to a stationary dark 
soliton, but the experiments are carried out in a harmonic trap [B^ and dissipative effects 



were shown to occur that accelerated the solitons 32,33 . 



C. Bright solitons 



For the case of attractive interactions (i.e. Qsc < 0), the kinetic energy of the BEG 
can be balanced by the nonlinearity yielding non-spreading wave packets. These solutions 
correspond to spatially localized bright solitons [ PU| , PU[ . The one-parameter solution to Eq. 
(O) then reads 



f{Z,t) 



n^/^/2-— sech 



n 




^i{4>{z,t,v)+^it/n) 



(16) 



with V the velocity parameter and /i, /q and Vq as in the dark soliton case. The phase 
given by 



IS 



0(Z,t,t;) 



k 



vt) ( V 



(17) 



Such ID bright solitons could in principle be realized in cigar shaped BEGs with negative 



scattering lengths, for example, in ^Li or ^^Rb by using Feshbach resonances p4|-p^ to 
tune the scattering length. To the best of our knowledge such experiments have not been 
attempted. In two or more dimensions negative scattering lengths can lead to catastrophic 
collapse in homogeneous systems for large enough particle numbers. However, in quasi- ID 



systems with strong transverse confinement the solitons are rendered stable [37 



D. Gap solitons in optical lattices 



So far in our discussion dark solitons arise for positive scattering lengths and bright 
solitons for negative scattering lengths. However, it is possible to extend the range of 



options, e.g. bright solitons with a positive scattering length, by considering a ID BEC in 



a periodic optical lattice |l38| , |39| . This situation was first discussed by Zhang et al. |^ for 
a scalar or single component condensate, and yields atomic gap solitons. Historically, gap 
solitons were first considered in nonlinear optics as arising from the combination of optical 
nonlinearity and a periodic spatial refractive- index distribution |§-0. As an introduction to 
our treatment of gap solitons in spinor condensates, we briefiy review the physics underlying 



gap solitons in scalar condensates, referring the reader to the literature for more details |^ . 
Consider the ID GPE (^) with an additional periodic potential, or optical lattice, 
produced by periodic optical fields and applied along the Z-axis, so that 

i^-K7 = ~^^^ + ^°Pi cos^{koptZ)(p + U\ip\'^ip . (18) 

Here Vopt determines the strength of the optical lattice, and kopt is the wave vector charac- 
teristic of the optical lattice. 

Let us first consider the linear eigensolutions of Eq. (|18]) neglecting the nonlinear term. 



Introducing the band-index n and the wave vector k within the band, we can write the 
atomic wave function in terms of the periodic Bloch functions Un,k{Z) according to Bloch's 



theorem [41 



<^(Z,t) = <^„,fc(Z)e-^^-'=*/^ = u„,fc(Z)e^('=^-^".'=*/^) . 

Using this expression in Eq. ([TsD with f/ = 0, the eigenvalue equation associated with the 
cosine potential then reads 

En,k'^n,k = -^ -7^ + Vopt COS^ {koptZ)iPn,k , (19) 



the eigenvectors and eigenvalues of which are well known [0. For our discussion the im- 
portant fact is the occurrence of an atomic band structure given by the dispersion relation 
En,k = f>'^n,k- Cousider an atomic wavepacket for a given band-index n and a narrow spread 
of wave vectors Ak around a central value fco that is slowly varying on the scale of the 
periodicity. Then we can write the atomic wave function as a product of the slowly varying 
envelope x{Z, t) with the fast oscillating Bloch part 

(^(Z,t) = x(^,t)¥..,fcoe-^^"■'=o*/^ (20) 

In this approximation the group velocity Vg and effective mass frieff for the wavepacket will 
be 

_ 1 dEn^k I 1 _ 1 ^ En.k I ,„^s 

in terms of which the approximate dynamics of the slowly varying envelope xi^, t) may be 
described by 

where we have re-introduced the nonlinearity. 



We take for illustration the case ko = 0, that is zone center, so that Vg = 0. Then Eq. 
(P^) becomes identical with Eq. ([T2|) , so that they have the same bright and dark soliton 
solutions but with m -^ ^^eff- The key idea is that the effective mass rrieff can assume both 
positive and negative signs by a suitable choice of the band index n and fco, whereas for 
the spatially homogeneous case the atomic mass m is always positive. So now, for example, 
bright solitons can arise for a positive scattering length if the effective mass is negative, 
the general condition for bright solitons being rUeffU < 0, and conversely for dark solitons 
rUeffU > 0. We remark that the applied optical lattice can modify the effective mass along 
the Z-axis, but the effective mass along the transverse dimensions is still m. This leaves 
open the question of higher-dimensional gap solitons, a point that we do not further address 



here [43 



The solitons formed in optical lattices are referred to as gap solitons 10-0. This name 
derives from the fact that the soliton energies actually lie in the energy gaps of the atomic 
band structure ^. The significance of gap solitons in the context of atomic BEC is that they 
provide a means to realize bright solitons, that is, spatially localized atomic packets, even 
for positive scattering lengths. This may have important applications, e.g. for the coherent 
transport of atoms. In the following section we further develop the theory of atomic gap 
solitons in spinor condensates to quantify some of their most important properties. 

IV. ATOMIC GAP SOLITONS 

We now consider a system of two condensates in different Zeeman sublevels coupled by a 
spatially periodic two-photon interaction. This spatially modulated coupling again leads to 
an atomic band structure for the spinor condensate, and hence by our general arguments of 
the previous section, to gap solitons. The advantage of multicomponent condensates is that 
magnetic fields can now be used to phase-imprint the two Zeeman components differently, 
and hence to manipulate and control the gap solitons to a high degree. We demonstrate this 
flexibility in Section V by simulating an atomic interferometer based on gap solitons. 

A. The physical model for a spinor condensate 

The two-component Bose-Einstein condensate interacts with two counterpropagating, 
focused Gaussian laser beams of equal frequencies ui but opposite circular polarizations. 
The optical dipole potential associated with the applied laser beams is assumed to provide 
tight transverse confinement for the BEC in the {X, Y) plane, thereby forming a cigar- 
shaped condensate of transverse cross-sectional area At- As in Section |T|, we confine our 
discussion to the one-dimensional dynamics of the BEC along the Z-axis for simplicity. 

In addition to supplying a transverse optical potential, the laser beams can drive two- 
photon transitions between different Zeeman sublevels of the atomic ground state. For 
illustrative purposes we consider the case of Sodium and the two-photon coupling of the 
Zeeman sublevels \ — 1) = \Fg = l,Mg = — 1) and |1) = \Fg = l,Mg = 1). For example, 
starting in the | — 1) state this process involves the absorption of a 0"+ photon from the 
right propagating laser beam followed by emission of a o"_ photon into the left propagating 
laser beam. We must of course assume that the excited states involved in the atom-field 



interaction are far-detuned from the applied laser frequency, a necessary requirement to 
avoid the detrimental effects of spontaneous emission. 

By restricting our attention to the coupled states | ± 1) the effective single-particle 



Hamiltonian for our model system can be written as [§J,^ 



H. 



p2 

£_z_ 

2m 



+ gh6' |1)(-1 



^2iKiZ j^ 



-1)(1| 



,-2iKiZ 



(23) 



where we have omitted constant light-shift terms. Here Pz is the atomic center-of-mass 
momentum operator along Z, Ki = uJi/c is the magnitude of the field wave vector along Z, 
g is & coupling constant between the ground and excited states characteristic of the atom 
and transition involved, 5' = V'^E'^ /h^5, with the detuning 5 = uji — Ua, £ is the laser field 
amplitude, and V is the reduced electric dipole moment for the 3S'i/2-3P3/2 transition. 

The second term of the effective Hamiltonian (p3D describes the effective coupling of the 
two Zeeman sublevels via the applied laser fields. The exponential terms exp(±2ii^;Z) arise 
from the fact that the two-photon transitions involve the absorption of a photon from one 
light field and reemission into the other. Introducing a spinor macroscopic condensate wave 
function ip{Z, t) = [fi{Z, t), ip^i{Z, t)]"^ normalized to the number of atoms N, and including 
the many-body effects via a mean-field nonlinearity, we obtain the coupled GPE equations 



ih 



dip 
'dt 



Heffip + U\ip\'^ip 



(24) 



where U = Anh asc/fn as in Section and [HI , 



m 



Iv^ip -|- Iv^-iP, and we have assumed 



that the magnitude of the self- and cross-nonlinearities are equal for simplicity. 

It is convenient to re-express Eq. (p3) in dimensionless form by introducing the scaled 



variables t = rtc, Z = zl^ and ipj = i^j^/p^ where 



gS" ^' 



tchKi 



Equations 



then become 



'dr \iJ- 



m 



-MV^ e^''''^ 



gfiS' 



U 



(25) 



-2ikiz 



+ sgnig6'/U)\ij\ 




(26) 



where M = gS'm/2hKf is a mass-related parameter such that ki = Kil^ = 1/2M. For a 
discussion of characteristic values for the case of a Sodium condensate we refer the reader 
toRef. ira. 



B. Soliton solutions 



The spatially modulated coupling between the optical fields and the condensate induces a 
single-particle band structure with regions of negative effective mass. As mentioned before, 
this leads to the possibility of bright atomic solitons even for repulsive interactions Ei 
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Approximate analytic expressions for these solitons can be obtained by expressing the spinor 
condensate components as 

^±i(^,r) = e±^'='^e-^-/^^V±i(^,r), (27) 

where the field envelopes 0±i are assumed to be slowly varying in space compared to l/h. 
Neglecting the second-order spatial derivatives yields the coupled partial differential equa- 
tions 




(28) 

where the choice ±1 = sgn{g5' /U). For a red-detuned laser and our choice of g this becomes 
±1 = sgn(f/). Neglecting the nonlinearity one finds a linear dispersion curve as sketched in 
Fig. ^. In this particular case it is the hnear coupling induced by the laser field that creates 
the avoided crossing at fc = and the negative effective mass. Aceves and Wabnitz [^ have 
shown that the full equations ( pH]) have the explicit two-parameter gap soliton solutions (see 
also Ref. §) 

sin(?7) ( e^^ + e^'^^\ ufr,^'^''l\ ±io- 



sech 6* =1= -— e 



"/?7V2 V e2^ + e±^'?y V 2 

P.HV)( ^'' + ^^''\\,,h(e±'-^)e^-. (29) 
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v/2 V e2^ + e±*''y V 2 



Here — l<f<lisa parameter which controls the soliton velocity, < r] < tt is a shape 
parameter, and 

^ = — 7sin(r7)(z — t>r) , cr = — 7cos(?7)(f2; — r) . (31) 

Since we are interested in creating bright solitons in the presence of repulsive interactions 
we restrict ourselves to sgn(f/) = -|-1, corresponding to the choice of the upper sign in the 
analytic solutions. The characteristic length scale associated with the solitons is Ic-, so that 
the approximate solitons ( p9| ) are valid for Kil^ = 1/2M >> 1. 



C. Soliton properties 

From the dependence of the hyperbolic-secant on ^ = —^sin{ri){z — vr) in Eqs. (p9D, 
we identify the gap soliton parameter v = Vg/Vji as the group velocity Vg of the soliton in 
units of the recoil velocity Vr = Ic/tc = hKi/m. Since — 1 < f < 1, the magnitude of the 
group velocity is bounded by the recoil velocity. From Eqs. (p9l), one can extract further 

9 



important soliton properties, such as the number A^^ of atoms in the sohton and the sohton 
width Ws = Wglc- Specifically, the number of atoms in a particular gap soliton is given by 

Ns = At I dZ[\MZ,t)\'' + |¥^-i(Z,t)|2] , (32) 

where At is the effective transverse area. 

We can gain further insight into the analytical structure of the gap solitons by taking the 



extreme limit 77 << 1 of Eqs. (p9|). Using the definitions (pT]) and returning to dimensional 
units, we have then 

^i(Z,0) = ^Msech{Z/Wo){-ire^^''+''''>^ , 

/37 V 2 

(/p_i(Z,0) = ^,/^sec/i(Z/Vro)(-l)"e*[(^-^')^+"l , (33) 

7 V 2 

where 

W, = 3.UW0 = 3.44 (k/lEZ\ , K = -^. (34) 

Here Wg = 3.44PVo is the soliton width, the factor 3.44 being the numerical conversion from 
the width of the hyperbolic secant to the 1/e^ width of the distribution, and K = k/lc a 
velocity-dependent wave vector shift. Obviously the width decreases with increasing 77 and 
V. The soliton atom number obtained by combining Eqs. (|3^ ) and (p3D, 

N, = 2{AtIcPc)v ■ (1 - v^) , (35) 

predicts that increasing the shape parameter 77 increases the atom number, and faster solitons 
have lower atom numbers. 

An essential point to keep in mind is that the gap solitons are coherent superpositions of 
the two Zeeman sublevels. The approximate solutions (^) contain important information 
on the phase and amplitude relations that need to be created between them to successfully 
excite and manipulate gap solitons. In particular, they show that there is always a spatially 
homogeneous tt phase difference between the two states. In addition, the two components 
have the spatial wave vectors 

K±i = K±Ki, (36) 

with K = —'-fv/lc- 

Finally, it follows from dividing the amplitudes of the two components that 



V9i(Z,t) 



(p-i{Z,t) 



2 



^-(^1- (^^' 



which shows that their relative occupation depends on the soliton velocity parameter v. 
For V = the sublevels are equally populated, but as f — > 1 the |1) sublevel has a larger 
population, and vice versa for i; — > — 1. 
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The characteristic time scale for the evolution of the gap solitons can be determined from 
the plane- wave exponential factors in Eqs. (p9|). Converting back to dimensional form the 
soliton period tg is defined as the time to accumulate a 27r phase, or in the limit rj ^ 

ts = 27rtcVl - v^ . (38) 

Physically, tg corresponds to the internal time scale for the gap soliton. In order to observe 
a soliton-like behavior, it is therefore necessary to investigate the atomic propagation over 
several periods. 

We conclude this section by noting that Aceves and Wabnitz [^ have shown that the 
gap soliton solutions (^9]) are stable solutions of Eqs. (^) in that they remain intact during 
propagation, even when perturbed away from the exact solutions. However, one should 
remember that Eq. (^81) is only an approximation to the exact system of equations (p6|), 
so that in general the gap solitons are solitary wave solutions only. As such, they are not 
guaranteed to be absolutely stable. 



V. MAGNETO-OPTICAL CONTROL 

Summarizing the main results of the preceding section, we have seen that gap solitons 
require the right population and density distribution in each Zeeman sublevel, a certain 
phase difference between these sublevels, and appropriate plane- wave factors e'^±i^. Based 
on these properties, we now show that a magneto-optical scheme involving a combination 
of pulsed coherent optical coupling and of phase-imprinting using spatially inhomogeneous 
magnetic fields turn out to be an adequate and convenient tool to manipulate these solitons. 

A. Manipulation Tools 

The shapes of the Hartree wave functions corresponding to the two Zeeman sublevels are 
hyperbolic-secant, which we approximate by a Gaussian in the following. They could for 
example be initialized in an optical dipole trap [^]. Manipulating the gap solitons reduces 
therefore to the problem of controlling the populations and phases throughout the spinor 
condensate. The coherent optical coupling can be achieved e.g. by a laser pulse of frequency 
uji propagating perpendicularly to the Z-axis and with linear polarization perpendicular to 
that axis. For sufficiently short pulses, one can neglect changes in the center-of-mass motion 
of the atoms during its duration, leading to a very simple description. We assume for sim- 
plicity a plane-wave rectangular pulse of duration tp and of spatial extent large compared 
to the soliton. The Hamiltonian describing the coupling between this pulse and the con- 
densate is then the same as in Eq. (|2BD, but without the linear momentum exchange terms 
exp{±2iKiZ) and the kinetic energy term, and with 6' —>■ 6p. The state of the system after 
the pulse is then easily found to be 

f Vii^p) \ = ( ^°s^ ^sinx\ / V5i(0) \ 
\V-i{tp) J \-ismx cosx / VV5-i(0) / 

.M.(x)(X;',y. (39) 
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where x = 9^'ptp is the excitation pulse area and the operator Ml can be used to control the 
population transfer by an appropriate choice of x- 

The required phase relationship between the two states can be achieved via Zeeman 
splitting. Considering for concreteness a spatially inhomogeneous rectangular magnetic field 
pulse of duration ts we have, neglecting again all other effects, 

ih-^±i{Z, t) = ±fLBgF {Bo + B'Z) ^±i{Z, t) ■ , (40) 

where gp is the Lande g-factor of the hyperfine ground state, hb is the Bohr magneton, 
Bq the spatially homogeneous component of the magnetic field, and B' its gradient, the 
direction of the magnetic field being along the Z-axis. The application of this field results 
in the state 

^i^+^BZ) 1 <^_,(0) ) 




^i(O) 
¥^-1(0) 



(41) 



where 'd = —{fiBgF/^)BotB and Kb = —{fiBgF/^)B'tB are the imprinted phase shift and 
phase gradient (or wave vector), respectively. That is, the application of the magnetic pulse 
results in a phase difference of 2-t} between the two Zeeman sublevels, and in addition it 
imparts them wave vectors -^Kb- 

B. Excitation and Application 

To illustrate how stationary and moving solitons can be excited using this scheme, we 
start from a scalar condensate in the | — 1) state, ^{Z, 0) = [0, (^o(^)]"^, with spatial mode 

with Ws and Ns the width and atom number of the gap soliton desired. This Gaussian is 
chosen to approximate the hyperbolic-secant structure of the analytic gap soliton solution 
in Eq. (|33|). For a stationary solution we need to prepare the Zeeman sublevels with equal 
populations and with a tt phase difference. We further need to impose wave vectors which 
are equal in magnitude but opposite in sign, K±i = ±Ki. This can be achieved by applying 
a laser pulse of area x = ^/4, followed by a magnetic pulse with {} = 7r/4 and Kb = Ki. 
The state then transforms as {t = tp + ts) 

ip{Z,t) = MB{T^/^,Ki)ML{T^/A)ip{Z,Q) 



g-T" / pi-KiZ 



^ I e ""e 



i-K „—iKiZ 



MZ) . (43) 



2 



Figure § shows the resulting stable evolution of the total density \'^{Z, t)\ . As a result of the 
Gaussian approximation to the exact solution there are some slow oscillations imposed on 
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the motion, but the solution remains centered at Z = and stationary over a time t = 200 
ms, much longer than the soliton period for the chosen values. We note that although 
Section |V A| treated the pulsed excitations in an impulsive manner to intuitively understand 



their action, our simulations do not make this approximation. The numerics confirm the 
accuracy of the impulsive approximation for the parameters at hand; a consequence of the 
fact that we consider pulse durations significantly shorter than the soliton period (|38|) . 

The excitation of a moving soliton is slightly more complicated since the velocity depen- 
dent wave vector K = —jv/lc in Eq. (^3]) is no longer zero and the two components have 
different populations, see Eq. (^7\). However, our simulations show that the proposed scheme 



performs well on this task, as well as on splitting solitons and reversing their directions |^ 

Solitons present some advantages for atom interferometry in that they are many-atom 
wavepackets which are immune to the effects of spreading, hence allowing longer path 
lengths, and also increased signal-to-noise for large atom numbers. Typically many-body 
effects limit the utility of high-density wavepackets due to spatially varying mean-field phase 
shifts, but solitons have the cardinal virtue that they have fixed spatial phase variations. 
Thus they may provide a key to making maximal use of high density sources for atom inter- 
ferometry. Indeed, due to these very properties they have long been advocated for all-optical 
switching applications. 

VI. SUMMARY 

In this paper, we have reviewed important aspects of the matter-wave solitons that can 
be launched in weakly interacting Bose-Einstein condensates. So far, only dark solitons 
have been demonstrated experimentally, but the rapid progress in the optical manipulation 
and control of condensates is expected to lead to the demonstration of bright gap solitons. 

As an illustration of the potential use of these objects, we consider a soliton-based non- 
linear atomic Mach-Zehnder interferometer. Our specific demonstration involves an initial 
scalar condensate that is split into two oppositely moving solitons along Z, see Fig. ^ for 
t < 60 ms. At t = 60 ms laser and magnetic pulses are applied which act to reverse the 
direction of the two solitons. This process causes some loss of atoms in both solitons, as 
before. The two reversed soliton components come together again at t = 120 ms. Since the 
colliding solitons are predominantly in opposite orthogonal Zeeman sublevels, the interfer- 
ence pattern appearing during the collision is due to the contamination of each soliton by 
the other state. Fig. |^ shows the interference at the soliton collision in the total density 
(solid line) and also the individual Zeeman sublevels, with a fringe contrast of around 30%. 
These results demonstrate the potential use of gap solitons for realizing nonlinear atom 
interferometers with high brightness sources. 
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FIGURES 
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FIG. 1. Density distribution (a) and spatial phase (b) of a stationary dark soliton with 6 = ir. 
The dip in the density has a width ~ ^o- The scheme for the generation of dark solitons by phase 
imprinting is shown in (c) 
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FIG. 2. Absorption images of BECs with dark soliton structures propagating along the lon^ 
condensate axis for different evolution times tp„ 
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FIG. 3. Linear dispersion curve (arbitrary units): Shown are the two branches of the dispersion 
relation exhibiting an avoided crossing at A; = (solid line). The negative curvature in the lower 
branch around k = defines a negative effective mass that enables gap soliton solutions. The free 
dispersion when no linear coupling is present corresponds to two displaced parabolas {dashed line) 
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FIG. 4. The propagation of the total density of a stationary sohton (v = 0) over around 50 
sohton periods. Peak density oscihations are due to imperfect initial conditions 



20 



-1000 




120 



Z[)im] 500 



80 
60 
40 t [ms] 



FIG. 5. Demonstration of an atomic Mach-Zehnder interferometer: The initial condensate is 
split into two counterpropagating solitons, then their direction is reversed and they collide (shown 
is the total density) 
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FIG. 6. Interference pattern when two solitons collide in the Mach-Zehnder configuration: 
The total density {solid line) is symmetric, whereas the interference pattern of each of the two 
orthogonal Zeeman states {dotted and dashed line) is asymmetric due to the shape of the colliding 
wavepackets. The contrast in the total density pattern is around 30% 
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